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In this paper, the density-matrix renormalization group method is employed to investigate the 
fractional quantum Hall effect at filling fractions v = 1/3 and 5/2. We first present benchmark 
results at both filling fractions for large system sizes to show the accuracy as well as the capability 
of the numerical algorithm. Furthermore, we show that by keeping a large number of basis states, 
one can also obtain accurate entanglement spectrum at v — 5/2 for large system with electron 
number up to N e = 34, much larger than systems previously studied. Based on a finite-size scaling 
analysis, we demonstrate that the entanglement gap defined by Li and Haldane p| is finite in the 
thermodynamic limit, which characterizes the topological order of the FQHE state. 

PACS numbers: 73.43. Lp, 73.43.Cd 
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I. INTRODUCTION 

In condensed matter physics, one of the major chal- 
lenges is to understand the strong correlation effect in in- 
teracting electron systems. The fractional quantum Hall 
effect (FQHE) systems^ are primary examples, where 
new quantum phases emerge with fractional quasipar- 
ticle excitations 3 - resulting from such effect. Theoretical 
understanding 3-7 of the odd denominator FQHE states 
in the lowest Landau level (n = 0) has been fully de- 
veloped soon after the experimental discovery 2 . The 
even denominator^ quantum Hall states in second Lan- 
dau level (n — 1) appear to demonstrate complex nature 
beyond the understanding of a unified theory, which are 
under intensive studies^—. 

Theoretical approache o 7 i 10 ~ 12 i 16 have predicted a va- 
riety of possible candidate states for even-denominator 
quantum Hall systems, some of which are exotic in na- 
ture with quasiparticles obeying non-Abelian statistics. 
However, very often such theories cannot provide solu- 
tions to microscopic theoretical models describing realis- 
tic electron systems. In this aspect, computational stud- 
ies have played an important role to determine the quan- 
tum state for such systems. Pioneer works have been 
done using exact diagonalization (ED) method to estab- 
lish the nearly perfect overlap of the Laughlin wave func- 
tion with the exact ground state wavefunction of small 
electron systems and provide microscopic understanding 
of the nature of the FQHE 4 — In recent years, ED has 
been widely used to study FQHE for pure systems 2 ^— 
as well as disordered systems^ to probe the nature of 
various quantum phases and transitions. It has been re- 
cently established that ED can also be used to identify 
the topological nature of a quantum phase based on en- 
tanglement entropy and the entanglement spectrum of 
the quantum systemsii^— . However, since the Hilbert 
space increases exponentially with the system size, ED 
is limited to systems with a small electron number (typi- 



cally restricted to electron number N e = 14 ~ 20 depend- 
ing on the Landau level filling number). This limitation 
becomes a severe problem when the finite-size effect is 
strong, which is usually the case for quantum states in 
the higher Landau levels or with even denominator filling 
fractions. 

The density-matrix renormalization group (DMRG) 
developed by White is a powerful method for studying 
interacting systems with accuracy controlled by the num- 
ber of states kept in DMRG blocks 3 ^ 1 . It has been 
widely applied to quasi-one-dimensional systems, provid- 
ing essentially exact results for spin or electron systems. 
However, its application to two-dimensional systems re- 
mains to be quite challenging. It is generally believed 
that the number of states desired to be kept in each block 
should grow exponentially with the increase of the sys- 
tem width to catch up the entanglement entropy between 
two coupled blocks 3 -. Nevertheless, FQHE systems can 
be modeled as one dimensional system with long-range 
Coulomb interaction, which may become accessible us- 
ing current computational power. Shibata and Yoshioka 
made the first attempt to develop DMRG algorithm to 
study quantum Hall systems in torus geometry 33 about 
ten years ago. Various fillings have been studied 3 -! by 
keeping hundreds of states obtaining useful information 
of larger systems at integer fillings or for compressible 
states. Recently, a great progress has been made by 
Feiguin et.a l 35 ^ 36 developing new DMRG algorithm keep- 
ing up to 5,000 states in their work. The hard problems of 
the incompressible FQHE at v = 1/3 and 5/2 have been 
extensively studied, providing convincing high accuracy 
ground state and excited state results for larger system 
sizes up to 20 and 26 electrons at v = 1/3 and 5/2, re- 
spectively. Moreover, DMRG has also been applied to 
study bosonic quantum Hall effect^ 7 .. 

In this paper, we study FQHE systems based on our 
newly developed DMRG code, which substantially im- 
proves the accuracy of the results through keeping more 
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states and managing DMRG process with higher effi- 
ciency. The benchmark results and error analysis are 
presented for v = 1/3 and 5/2. For the same system 
sizes studied before, we reduce the error substantially es- 
pecially for 5/2 systems, while we access systems with 
more electrons with high accuracy. We also obtain ac- 
curate entanglement spectra for v = 5/2 systems, which 
identify the topological order at larger system sizes. The 
entanglement gap for v — 5/2 FQHE first revealed by Li 
and HaldaneA based on ED calculation, remains finite in 
the thermodynamic limit established through finite-size 
scaling analysis. 

II. HAMILTONIAN AND METHOD 

In this paper, we study the FQHE in the spherical 
geometry^, where electrons are confined on the surface of 
a sphere with radius R. The total magnetic flux through 
the spherical surface AnR 2 B are quantized to be an inte- 
ger 2S multiple of the flux quanta. Assuming that elec- 
trons are polarized by the magnetic field and neglecting 
the Landau level mixing, the Hamiltonian in the spheri- 
cal geometry can be written as 

U = ^ X! ( f nim2\V\m 3 m 4 )al ril al ri2 a m3 a m4 (1) 

mi ,7712 ,7773 > m 4 

where the rrii is the z-component of the angular momen- 
tum and rrii — —L, —L + 1, • • ■ , L, with L = S + n being 
the total angular momentum and n is the Landau level 
index. V is the Coulomb interaction between electrons 
in units of j^, with = y/hc/eB being the magnetic 
length. In the next sections, we also use the renormal- 
ized magnetic length 1^°°^ to rescale the energy unit if it 
is mentioned explicitl y 22 ! 24 . a m (a] T1 ) is the annihilation 
(creation) operator at the orbital m. The problem can 
be naturally studied using the DMRG method in mo- 
mentum spaced as nonzero matrix elements in the Eq. 
([1]) only exist between orbitals satisfying the angular mo- 
mentum conservation relation mi + m2 = 7713 + TO4. 

In our DMRG process, we arrange the 2L + 1 or- 
bitals into a one-dimensional chain, corresponding to 
states with z-component of the angular momentum m 
taking different values L,L— 1, • • • , — L + l,— L from left 
to right. In the initial process, we start from a small 
system with two blocks (each block has only one or- 
bital) and two single orbitals in the middle forming a 
configuration Bl • •Br- We write the Hamiltonian ([1]) 
terms using the product of block operators and single 
site operator o 30 ' 31 ' 38 and diagonalize the Hamiltonian to 
obtain the ground state. The reduced density matrix 
is formed and diagonalized by following the standard 
DMRG22 procedure. A new block Bl (Br) is formed 
by including the single site into the current block. This 
procedure is repeated until the ID chain grows into the 
desired system size with 2L+1 orbitals. From this point, 
the finite-size sweeping algorithm is employed until we 
obtain the converged results. 
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FIG. 1: (color online) Error in the ground state energy per 
electron 5eo is shown at (a) v = 1/3 and N e = 20, with 
1,600 (O), 3,200 (□), and 5,000 (0) states kept, respectively. 
The triangle data(obtained by keeping 5,000 states) illustrate 
a situation where DRMG runs into a local minimum during 
initial 12 sweeps, see text for details. The reference ground 
state energy is -0.4210509, obtained by DMRG with at most 
24,000 states kept, (b) v = 5/2 and N e = 24 with 1,600(0), 
3, 200(D), and 5,000(0) states kept , respectively. The data 
showed by x are obtained by keeping 8,000, 10,000 and 14,000 
states for each 4 sweeps starting from the first sweep, and 
16,000 states for the rest of sweeps. The reference ground 
state energy is -0.3876150, obtained by DMRG with at most 
30,000 states kept. 



Using our momentum space DMRG, the Hamiltonian 
([I]) is diagonalized in the sector with fixed electron num- 
ber N e and the total angular momentum z-componcnt 
L z . In the initial process of the calculation, only a small 
fraction of the momentum sectors in each block con- 
tributes significantly to the ground state wavefunction. 
States in other sectors, which have zero eigenvalues in 
reduced density matrix, are being discarded. However, 
these discarded states may become also important to the 
ground state wavefunction in the later stage when system 
grows to the full length 38 . To overcome this problem, we 
need to keep additional sectors from very beginning. For 
this purpose, we set a minimum number of sectors iV™ n , 
which in practice is about 3^5 times larger than the 
number of sectors after convergence. In each of selected 
sectors, we keep at least two states with largest eigen- 
values in the sector. The remaining states are selected 
following the standard DMRG procedure to minimize the 
truncation error. In the sweep process, N™ n is gradually 
decreased as a function of sweep number. 

In Fig. [Ha), we show the error of the obtained ground 
state energy per electron at v ~ 1/3 for different number 
of states kept in each block with respect to the fully con- 
verged reference energy obtained by keeping much more 
states. A reasonable accuracy around 10~ 4 is reached 
with the electron number N e = 20 (bigger than the 
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largest ED size by six electrons) by keeping only 1,600 
states, while we achieve the accuracy of 6 x 10~ 6 by 
keeping 5, 000 states. However, the convergence becomes 
more difficult for v — 5/2 system, as we show in Fig. 
QJb). For N e — 24 (bigger than the largest ED size by 
four electrons), with m = 5, 000 states kept, the accuracy 
we can achieve is about 2 x 10 -4 , comparable with the re- 
sults obtained by Feiguin et al2£. Further increasing the 
number of states kept to 16, 000, we are able to reduce 
the error by a factor of five. 

Controlling N™™ also allows us to overcome the lo- 
cal minimum trapping, which could trap the DMRG ob- 
tained state in an excited state. As we show in Fig. Ufa) 
by triangles, the energy is pinned to a local minimum 
from the fifth sweep to twelfth sweeps, and the error is 
much larger than the data showed by <0 starting from a 
different initial state. To overcome the problem, we sim- 
ply increase N™ n at the twelfth sweep to allow a larger 
number of sectors to get into the Hilbert space. These ad- 
ditional sectors bring in significant quantum fluctuations 
and eventually get the state out of the local minimum. 



III. RESULTS 

A. Ground state energies at v = 1/3 and 5/2 

To demonstrate the accuracy of our DMRG calcula- 
tions, we have obtained the ground state energy for sys- 
tem up to 24 electrons at v = 1/3 by keeping up to 
20, 000 states, which leads to a truncation error smaller 
than 10 -11 in the final sweep. The maximum dimen- 
sion of the Hilbert space diagonalized is of the order of 
2 x 10 7 . In Fig. [21 we show the ground state energy 
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FIG. 2: (color online) (a) Finite-size scaling of ground state 
energy per electron at v = 1/3 with a quadratic extrapolation, 
(b) Finite-size scaling of ground state energy per electron at 
v = 1/3 with a linear extrapolation. Energy unit is renormal- 
ized by . The error bars in both figures are much smaller 
than the size of the symbols. 



as a function of l/N e with 2L = 3(N e — So), where a 
shift of So = 1 has been used. In the upper panel nu- 
merical data are extrapolated to thermodynamic limit 
by a quadratic function of 1/-/V e , which gives the ground 
state energy per electron eo = —0.41016(2). In the 
lower panel we rescale the energy by the renormalized 
magnetic lengt h 22 ' 24 and extrapolate the numerical data 
linearly, leading to e = —0.410145(15) demonstrating 
consistency between two extrapolating methods. While 
our results are essentially in agreement with previous re- 
sults recently obtained by Feiguin et al2£ using DMRG 
method, our accuracy is improved by keeping much more 
states, which opens opportunity for studying larger sys- 
tems. Now we turn to the study of the FQHE at v = 5/2. 
We calculate the ground state energy up to N e = 34 elec- 
trons, with at most 24,000 states kept. The maximum 
dimension of the space we diagonalize is around 3.5 x 10 7 . 
The truncation error is of the order of 10~ 7 indicating 
larger error in DMRG comparing to v = 1/3 case due to 
larger entanglment entropy in the ground state as well as 
smaller gap separating the excited states from the ground 
state. The data at the largest system size were obtained 
within three weeks on one 12 cores Xeon server. In upper 
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FIG. 3: (color online) (a) Finite-size scaling of the ground 
state energy at v = 5/2 with a quadratic fitting, (b) Finite- 
size scaling of the ground state energy at v = 5/2 with a 
linear fitting. Energy unit is renormalized by Iq and only 
N e > 18 is used in view of the strong finite-size oscillation. 

panel of Fig. [31 we show the ground state energy e as 
a function of 1/N e . With a quadratic fitting, the extrap- 
olated value in the thermodynamic limit is —0.3645(6). 
In the lower panel, we rescale the ground state energy by 
the renormalized magnetic length /q 00 "* and extrapolated 
ground state energy is —0.36415(45) by a linear func- 
tion. We further calculate the excitation gaps at v = 5/2 
up to N e = 26, including neutral exciton ga.p(A exc ) and 
the charged excitation gap(A), following the definitions 
in Ref. [24j. Some of the data are present in Table [J as a 
function of N e , where the excited states "aliased" to other 
quantum Hall states are exclude d 23 ' 24 . The estimated 
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TABLE I: Excitation gaps at v — 5/2 as a function of N e 
Some data for exciton gaps are not shown here. 



the following, we demonstrate the success of our DMRG 
calculations in this aspect by keeping up to 24,000 states. 
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FIG. 4: The low-lying entanglement spectrum of N e 
the partition P[0|0] obtained by DMRG are shown. 
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gaps in the thermodynamic limit are A exc = 0.032±0.004 
and A = 0.029 ± 0.003, in agreement with early result, 
although we observe that the strong finite-size oscillation 
persists even for our larger system size data. 



B. Entanglement spectrum and entanglement gaps 

By dividing 2L + 1 orbitals into two parts, A and B, 
the ground state can be written, according to Schmidt 
decomposition, as 



- £4/2 lv4)®Ki), 



(2) 



where the singular eigenvalues exp(— &/2) obtained from 
diagonalizing the reduced density matrix define the en- 
tanglement spectrum £j. The \ip l A ) and \ij) l B ) are the 
orthogonal basis states of part A and B, respectively. 
In the pioneer work of Li and Haldanei, the entangle- 
ment spectrum for system at v = 5/2 has been ana- 
lyzed based on exact diagonalization calculation up to 
N e = 16 electrons. The obtained entanglement spec- 
trum show the same structure as the conformal field the- 
ory (CFT) for the Moore-Read state, below an entan- 
glement gap while the non-CFT type of spectrum ex- 
ist above the gap. Therefore, the entanglement spec- 
trum reveals more information than the entanglement 
entrop y 26 i 27 . Although the entanglement spectrum and 
entropy are naturally obtained in DMRG, they are much 
harder to converge^ than the ground state energy. In 
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FIG. 5: The low-lying entanglement spectra of N e = 34 are 
shown for three partitions P[0|0], P[0|1] and P[l|l]. 

We note that our two parts of system A and B used 
for Eq. (2) include one block and one site defined for the 
DMRG process, respectively, thus, the quantum numbers 
from both parts satisfy the equations N^+N® = N e and 
Lf + Lf = (the ground state for 5/2 FQHE has total 
L z = 0). The entanglement spectrum can be labeled by 
and Lf in part A. We recall thati^ for the ground 
state at v = 5/2, the highest-density "MR root configura- 
tion" has a pattern of "11001100. . .110011" correspond- 
ing to the "generalized Pauli principle" that no group of 
4 consecutive orbits contains more that 2 particles. Con- 
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sequently, there are three distinct ways of partitioning 
the orbits as between two zeros, between zero and one, 
and between two ones. These three different partitions 
of the root configuration, identified as P[0|0], P[0|1], and 
P[l|l] can be easily obtained in the sweeping process^ 
during the DMRG calculations. 

To check the accuracy we compare the entanglement 
spectrum obtained by DMRG with that obtained by ED 
for N e = 16 for the partition P[0|0], as shown in Fig. [4] 
DMRG results precisely reproduce ED results 1 . Next, we 
turn to consider larger size. In Fig. [51 we show the en- 
tanglement spectra of three partitions for N e = 34. The 
low-lying spectra have the identical counting structure 
viewing as a function of the AL = Lf max — Lf (here the 
L^max — 288 for the partition P[0|0]) comparing to the 
results presented by Li and Haldane [TJ], establishing the 
identical topological order as in the Moore- Read state at 
this larger system size. Furthermore, the same spectrum 
structure has also been obtained for a wide range of even 
electron numbers 8 < N e < 34. 
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FIG. 6: (color online). Entanglement gap at AL = vs 1/N e 
for three different partitions. Dashed blue line is a linear fit 
for iV e > 16. 

Finite-size calculation shows that the CFT part of en- 
tanglement spectrum is protected by an entanglement 
gap from the non-CFT part and, the entanglement gap 
plays the key role of measuring the robustness of topo- 
logically ordered stated. In literature, the entanglement 
gap is obtained by ED and limited to a small sizei^ 
and whether it persists in the thermodynamic limit re- 



mains an open question. Thus, identifying whether the 
entanglement gap persist in the thermodynamic limit is 
particularly important for the establishment of the topo- 
logical^ ordered state in such systems. This is different 
from the finite-size wavefunction overlaps between the 
wavefunction of a realistic interacting system and model 
wave-functions (Laughlin or Moore-Read wavefunctions) 
as such overlaps have to go to zero with the increase 
of the system size related to the presence of the generic 
entanglement spectrum in the system. Here we have ob- 
tained the entanglement spectrum gap for a relatively 
large range of the system sizes, which allows us to ex- 
trapolate the gap to the thermodynamic limit. In Fig. [6j 
we show the entanglement gap 8q at AL = as a function 
of 1/N e for three partitions. Entanglement gap oscillate 
strongly for small sizes, and the oscillation magnitude 
becomes smaller as N e increases. Moreover, the lines for 
the three partitions almost fall into one curve. We fit 
the entanglement gap as a linear function of 1/N e and 
neglect higher order corrections due to the oscillation of 
the data and limitation of the size. It yields So = 4.6±0.6 
in the thermodynamic limit and thus our work confirms 
the conjecture of a finite entanglement gap made by Li 
and Haldanei. 



IV. CONCLUSIONS 

We have presented systematic numerical results ob- 
tained by a newly developed DMRG program for FQHE 
systems. We have substantially improved the DMRG al- 
gorithm and obtained accurate results for ground state 
energy and excitation gap for larger systems than previ- 
ous works by ED and DMRG at v — 1/3 and 5/2. In par- 
ticular, we demonstrate the robustness of the low-lying 
CFT entanglement spectrum for large system size with 
N e = 34 and the finite entanglement gap in the thermo- 
dynamic limit at v = 5/2 based on finite-size scaling. 
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